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ABSTRACT 

We present calculations for a magnetised hybrid wind model for Betelegeuse 
(a— Orionis). The model is a direct application of our previously derived theory, com- 
bining a canonical Weber-Davis (WD) stella r wind with dust grains in the envelope of 



an AGB star (see Thirumalai & Heyl 



2010). The resulting hybrid picture provides a 



mec hanism for solving the problem of lifting stellar m aterial up from the photosphere 



(e.g. Harper et al.|2009 Guandalini et al.]|2006 Jura [1984) and into the circumstellar 
envelope. It also predicts wind velocities in agreement with current estimates. Our 



approach r eveals that magnetic fields in supergiant stars like Betelgeuse (seejAuriere 
et al.|2010| ), may play a vital role in determining the nature of the stellar outflow and 
consequently, opens a new avenue of investigation in the field of hybrid stellar winds. 

Key words: Betelgeuse, a— Orionis, AGB wind, Weber-Davis wind, MHD, dust- 
driven wind. 



1 INTRODUCTION 

Betelgeuse (a— Orionis), one of the closest cool-evolved su- 
pergiant stars to Earth, has been the focus of research 
for many decades now. However, despite the considerable 
amount of attention it has received, our understanding of the 
wind of the star and the nature of its mass efflux remain as 
yet, mysterious in more ways than one. Even though modern 
understanding of the stellar wind of this M-type supergiant 
star, despite formidable challenges on both the observational 
and theoretical fronts, has progressed well beyond the rudi- 
mentary stages (see Holzer fc MacGregorj 1985;, David fc Pa- 
poular|[r9 90 Dan chi et al.||1994| [Smith et al 



Brown||2006i |Harper et al.| 2008) Plez fc Lambert 20021 



2oom 



Harper 



Bernat et al.||1979| |Noriega-Crespo et al.||1997| |Ueta et al.| 
2003 1, still however, significant gaps in our understanding re- 



main quite entrenched. For example, observations in neither 
the ultra-violet (UV) nor the infra-red (IR), reveal blue- 
shifted emission lines; a requisite signature of gas outflow 
(see Harper et al.|2069 1 , should the wind get sufficiently ac- 
celerated close to the star. It is also still not known how 
the stellar wind is supported. We do however understand 
that dust forms at small and large distances from the pho- 
tosphere, at radii r ^ 1.33-Ro (e.g. Perrin et al. 20071 and 
r > 25 - 30-Ro (e.g. [Bester et al.||1991| |1996[ [Danchi et al 



1994 Harper et al. 2009 Skinner et al. 19971 respectively 



where Rg is the optical photospheric radius. While there is 
now evidence of dust forming species such as alumina being 
present in small quantities, in the lower reaches of Betel- 



geuse's atmosphere (~ 1.3 — 1.5Ro) (see Perrin et al.||2007 



E- maibanandQphas.ubc.ca (AT); heyl@phas.ubc.ca (JSH) 



Verhoelst et al.|2006 1 , there is no evidence that this results 
in a dust-driven wind at these distances (e.g. Harper et al. 
2009 1. It has been conjectured that alumina may form close 



to the star and then, once transported to large distances 
(r Z 30Ro), may provide nucleation sites for silicate dust 
to form (see Verhoelst et al.|[2006 l. However, this scenario 
still requires transport of stellar material to these larger dis- 
tances in the atmosphere and presently, it is not known how 



this is achieved (see Harper et al. 2009 Guandalini et al 



2006[|Jura|,1984) . It is believed that a combination of MHD 



effects, pulsation and convection may be responsible, but 
there are presently no models that demonstrate this un- 
equivocally (see [Harper et aL 20091. On the other hand. 



efforts at modelling the stellar wind of Betelgeuse (and red 
supergiants and AGB stars) remain almost as disparate fac- 
tions, modelling either solely magnetohydrodynamic (MHD) 
or acoustic waves (e.g. 'Hartmann & MacGregor'1980| 



Hart- 



mann fc Avrett 1984 Falceta-Gongalves & Jatcnco-Pereira 
2002 Falceta-Gongalves et al.|2006 l or relying upon a com- 



bination of dust- (see Elitzur fc Ivezic|2001 Liberatore et al. 
2001 Woitke fc Niccolini|2005 1 and pulsation-driven mecha- 



nisms (see [Elitzur et al.|2 003':'Bowen'1988':'Bowe n fc Willson 
T99T| [Fleischer et al.|1992 . Wood 1979 Bcdijn 1988[[Gail fc 



Sedlmayr[[1988 |. Each of these models achieves a modicum 
of success, while ignoring stellar rotation altogether. Though 
it is widely acknowledged, that both magneto-centrifugal ef- 
fects and the presence of dust grains in the gas, may be 



part of a greater coupled picture (see 


Thirumalai fc Heyl 


2010[[Guandahm et al.|2006[ 


Falceta-Gongalves fc Jatenco- 


Pereira 2002[ Falceta-Gonga 


Ives et al. [2006 Schroder fc 


Cuntz 


2007 


1, there have however been no efforts, thus far. 



at combining the two for a red supergiant like Betelgeuse. 
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In addition to these concerns, it is also clearly seen 
that the atmosphere around Betelgeuse is not spherically 
symmetric. Indeed, there are some inhomogeneities that are 



seen, attributable to clumps in the outflow (see Plez & Lam- 
bert|2002[ ). These concerns are compounded by the fact that 
observations reveal that the temperature structure has a 
complicated form. Differences are seen at the same radius 



in different parts of the atmosphere (see Lim et al. 19981, 



revealing a slight departure from a spherically symmetric 
picture. With regard to seeing a clear signature of a stellar 
wind in the form of blue shifted emission lines, the observa- 
tions have had severe difficulties in being able to resolve the 



spectra at close distances from the photosphere (see Harper 
et al.|2009 1. Indeed, observations have revealed that in parts 
of Betelgeuse's atmosphere there is redshifted emission indi- 
cating that there may be infall of matter back onto the star 
(see |Lobel fc Dupree|2001[ ). 

The rotation rate of Betelgeuse has also been reported 
in the literature (e.g. [Uitenbroek et aL]|1998| ). More re- 
cently, Zeeman observations also reveal unambiguously, that 
there exists a magnetic field on the surface of Betelgeuse 
(see |Auriere et al.|[20Tol [Petit et al.||2011[ ). These observa- 
tions naturally raise several questions regarding the role that 
magneto-rotational effects might play in the star's wind. 
Thus overall, we see that there are several unanswered ques- 
tions regarding the nature of the outflow from Betelgeuse. 
Here, we present the very first model integrating MHD and 
rotational effects with a dust-driven wind scenario for Betel- 
geuse. The aim of the current study is to shed a new light on 
these issues and attempt to answer, at least some of these 
concerns. In this study, we find that the presence of a small 
magnetic field, on the order of what was recently discov- 
ered; of about 1 G (see Auriere et al.|[2010| ), is sufficient to 
drive material from close to the stellar surface up and out 
of its gravity well, by means of a magneto-rotational wind. 
Dust condensation later occurs at large distances from the 
star (~ 30-Ro) which finally results in a hybrid-MHD-dust- 
driven outflow. This model provides a possible alternative 



resolution of this issue (see Harper et al. 2009 Guandalini 



et al. 2006 Jura 1984 1. It is to be mentioned at the very 



outset, that this model is merely suggested as an additional 
mechanism that can play a role in transport of stellar mate- 
rial and is not intended to supplant the altogether feasible 
models involving MHD or acoustic waves as likely candi- 
dates. 

The theoretical model that we developed for intermedi- 



ate mass Eisymptotic-giant-branch (AGB) stars (see Thiru- 



malai & Heyl 20101 has been extended to tackle the case 
of a supergiant like Betelgeuse. The reader is referred to 
our earlier work (see Thirumalai fc Heyl|2010 1 for details of 
the model. Here we shall present only the salient features 
of the theory. We also find that the wind velocities that are 
obtained from the model are in good agreement with cur- 
rent estimates (e.g. Harper et al. 2009 ; Falceta-Gongalves & 
|Jatenco-Pereira|2002| |Plez fc Lambert||2002f ! 



2 THE HYBRID- WIND MODEL 

For the purposes of modelling, there are some marked differ- 
ences that Betelgeuse exhibits from a run-of-the-mill AGB 
star. For example, Betelgeuse is much more massive (~ 



15Mq, Smith et al.||2009 l than an intermediate mass AGB 
star. It also has a far more extended and cooler atmosphere. 
More importantly, the dust condensation radius is much far- 
ther out, in terms of stellar radii, from the photosphere in 
comparison to a canonical AGB star. In the former, the 
primary process governing mass-loss is the radiation pres- 
sure on dust grains that form at close distances from the 
photosphere, typically no greater than 10i?o, coupled with 
strong stellar pulsation. The presence of a magnetic field 
then has the effect of playing a secondary role governing the 
gas dynamics. In an AGB star, the hybrid-MHD-dust-driven 
mechanism is such, that without the onset of dust forma- 
tion, it is not possible to achieve an outfiow and thereafter, 
the predominant energy exchange in the wind is between the 
magneto-rotational and gravitational components (see Fig. 5 
of [Thirumalai fc Heyl|2"oTo| . In such a scenario, it was seen 
quite crucially, that dust formation must occur prior to the 
hybrid model's sonic point and concomitantly, it was seen 
that the dust parameter, F^j was required to be less than 
unity for achieving a hybrid wind. Such is not the case for 
Betelgeuse where observations find only a small amount of 
dust inside r ;S 25 — 30 iZo- 

In our earlier study, we had additionally investigated 
the possibility of locating the dust condensation radius out- 
side the fast Alfven point. It was shown that such a hybrid 
wind is entirely theoretically possible, given typical parame- 
ters of an AGB star (see Fig. 9 of [Thirumalai fc Heyl|20"l0l ). 
It was concluded therein that while such a scenario is un- 
likely for a typical AGB star it may well apply to an alto- 
gether different type of stellar wind. It is this second type 
of hybrid-MHD-dust-driven model that is adapted herein to 
formulate a stellar efflux scenario for Betelgeuse, with dust 
formation occurring at a large distance ~ 30 Ro, from the 
star. 

Our theoretical model can be summed up as follows (the 



interested reader is referred to Thirumalai fc Heyl ( 2010 1 for 



the steps involved in the derivation). We imagine that we are 
looking down upon the two-dimensional equatorial plane of 
Betelgeuse. Therein, we assume that the magnetic field and 
the gas (and dust) velocity are functions of purely the ra- 
dial distance from the centre of the star. The poloidal (co- 
latitudinal) components of these vectors vanish; this is the 
fundamental assumption behind the canonical Weber-Davis 
(see Weber fc Davis|1967 hereafter WD) model for our sun. 
The gas forms the first fluid and carries the magnetic field. 
However, unlike the sun, in the atmosphere of Betelgeuse, 
embedded within the gas is a second fluid; the dust. The 
two fluids co-exist and are coupled to each other through 
drag. In an evolved supergiant like Betelgeuse, the circum- 
stellar atmosphere is cool enough (;£ lOOOK) that dust grains 
can condense out of the surrounding gas (e.g. Bester et al. 
T99T[ [T996| [Danchi et al.||1994[ [Harper et al.||2009| [Skin " 



ner et al. 19971. In fact, direct imaging has enabled esti- 
mates for the inner dust shell temperature to be ~ 700 K at 



around ~ 30-Ro (e.g. Bloemhof et al.|1984 |. Stellar radiation 
from the interior impinging upon the dust grains, can im- 
part enough momentum to power these exiguous solar sails 
and propel them outward through the surrounding gaseous 
matter. However, whilst moving through the gas, the dust 
drags the gas along with it, resulting in a prodigious and 
combined mass efflux of both dust and gas, from the star; 
the so-called dust-driven wind. It is to be kept in mind, that 



© 2011 RAS, MNRAS 000,[T|fT3| 



Betelgeuse Hybrid Wind 3 



the dust-to-gas mass-ratio is small; for Betelgeuse, it is ex- 
pected to be on the order of « 6 x 10~* — 5 x 10~^ or so 
(see Harper et al.|20dT |. Additionally, the dust grains in our 
model are assumed to be spherical in shape, thus presenting 
a circular cross-section for radiation pressure to act upon. 
The surface temperature of Betelgeuse (To) plays a key role 
in determining the hybrid wind parameters, such as the bulk 
gas velocity at the base of the wind and the location of the 
critical points (see below). We also assume that the gas has 
the thermodynamic equation of state of a polytrope, with 
a polytropic index 7 > 1, where a value of unity represents 
the isothermal limit. We model the fluid flow as an inviscid 
one and the electrodynamic properties of the fluid are taken 
to obey ideal MHD, i.e., there is no Lorentz force acting on 
the fluid and the electric and magnetic forces balance each 
other completely. Additionally, the model requires as input, 
a number of observed quantities, such as the mass of the star, 
its rotation rate, the surface magnetic field strength and the 
mass-loss rate. These and other parameters for Betelgeuse 
are listed in Table 1. With these ingredients, upon examin- 
ing the Euler equations for fluid flow for both the dust and 
the gas that are coupled with each other and invoking mass 
and energy flux continuity and ensuring that the divergence 
of the magnetic field explicitly vanishes within the governing 
equations, we can arrive at a steady-state description of the 
hybrid dual-fluid wind in the equatorial plane of Betelgeuse 
(see [Thirumalai fc Heyl||2010[ ). The radial equation for the 
gas velocity profile is then given by, 

dw w N{w, x) 
dx X D(w,x) 

where, w = u/ua is the gas speed normalised using the 
Alfven speed and x — r/vA, is the radial distance expressed 
in units of the Alfven radius. Hereafter, the subscript 'A' 
refers to values of the different variables at the Alfven radius. 
The quantities N{w,x) and D[w,x) are the numerator and 
denominator respectively and are given by. 



Table 1. Various parameters for modelling Betelgeuse. 



(1) 



N(w,x) = (^2'ySTi 
x{wx'^ - 1)^ -f- Snx^{w 
and 



-ii~rd-e{x~xa)) 



1) (1 - 3wx^ + {wx'^ + l)w) (2) 



D{w, x) = (10^ - jSriwx'^y-^') (wx'^ - 1)-^ - Snx^ x 



(3) 



In the above equations, the parameters St 



2kTA 
5- 



Sg 



; and Sn — — 5^ along with 7 uniquely determine the 

locations of the critical points, and hence the morphology of 
the family of solutions of Eq. ([TJ . Here Ta is the gas temper- 
ature at the Alfven radius, k is the Boltzmann constant and 
nip is the mass of a proton. The critical points are, as usual, 
defined as the locations at which both the numerator and 
the denominator vanish, thereby keeping the right-hand side 
of Eq. ([T| finite, these are the sonic point, the radial Alfven 
point and the fast point (e.g. Weber & Davis 19671. The 
presence of the Heaviside function in Eq. (|2|) represents the 
formation of dust at the location x = Xd, the dust conden- 
sation radius in units of the Alfven radius. The critical wind 
solution of Eq. ([l]) will yield the gas velocity profile, thereby 
enabling the determination of all other dependent variables. 



Parameter 


Symbol 


Value / Comment 


Mass 


M, 


~ 15M0 


Radius 


Ro 


~ 65OR0 


Mass loss rate 


M 


~ 3 X lO-^MQ/yr 


Surface magnetic 






field strength 


Bo 


~ 1 G 


Bulk surface gas 






velocity (radial) 


UO 


~ 10 fosc.O 






(vanishingly small) 


Surface temperature (effective) 


To 


~ 3650K 


Stellar rotation rate 


n 


~ 1.2 X 10"® rad/s 


Surface escape velocity 




9.39 X lO'^ cm/s 


Polytropic exponent 


7 


> 1 


Alfven Radius 


rA 


~ 25Ro 


Alfven speed 


UA 


~ 0.15t)osc,0 


Dust Parameter 


Td 


varied 


Dust grain radius 


a 


spherical grains 


Dust grain mass 


md 


~ A/Sna^Pd 



such as the dust velocity profile (to be discussed below), the 
Mach number as a function of distance from the star, the 
azimuthal velocity of the gas, the azimuthal component of 
the magnetic field, the temperature profile and the density 
structure of the gas in the envelope of the Betelgeuse. The 
critical solution of Eq. ([T]) is defined as one that starts off 
at the base of the wind sub-sonic, passes through the three 
critical points in a continuous manner and emerges super- 
Alfvenic at large distances from the star. The dust velocity 
profile is then given by (see [Thirumalai fc Heyl|2010[ ), 

/ I 2 \ 



v{r) = u{r) + 



a-th 



\ 



(4) 



where ath is the thermal speed given by ath = \j2kT j [im^ 
and /irriu is the mean molecular mass of the gas and is the 
dust grain number density, which is assumed to be given by, 
ndmd/ p ~ (S) where {5} is the average dust-to-gas ratio in 
the wind (see |Thirumalai fc Heyl|2010[ |. The dust in the cur- 
rent theory is treated in a simplistic and idealised manner, 
without rigorously considering the effects of dust radiative 
properties or including the effects of scattering and absorp- 
tion on the radiation pressure mean efficiency. While such an 
analysis would no doubt portray a more complete picture, 
the current rudimentary treatment nevertheless captures the 
salient features of the coupled outflow from the star. As our 
purpose here is to illustrate the feasibility of a hybrid-MHD- 
dust-driven wind model for Betelgeuse, the current simplis- 
tic treatment of dust was considered sufficient. 

Eq. |l]) is solved numerically, and the reader is referred 
to our earlier work (see [Thirumalai fc Heyl||2010| ) for com- 
plete details on the numerical methodology. The pertinent 
points of the method are conveyed below, in brief. 



3 NUMERICAL METHOD 

Eq. ([TJ is solved using the package ODEPACK employing 
a finite difference method with chord iteration with the 
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Jacobian supplied (see Hindmarsh 1983 Brown & Hind 
marsh|[l989 1. Initial conditions were supplied at the begin- 
ning of the integration. Typical error tolerances for con- 
vergence testing that were employed were on the or der of 
for bo th the absolute and relative errors (see 



10 



marsh 



Hind- 



19831). For a typical integration, step sizes of 10~^ 
or 10~^", in units of the Alfven radius, were employed de- 
pending upon the region of integration being near the crit- 
ical points or sufficiently away from them. This resulted 
in typically lO'' — 10^^° function evaluations. In the cur- 
rent study, in contrast to Thirumalai & Heyl (20101, we 



located the radial Alfven point at around 25i?o with an 
Alfvenic temperature of Ta ~ 720 K, and Alfvenic veloc- 
ity UA ~ 0.15wosc,o ~ 14 km/s, because we wanted to have a 
dust formation temperature of ~ 700 K with dust conden- 
sation occurring at ~ 30-Ro- The polytropic exponent 7 was 
varied and the locations of the sonic point and the fast point 
were found according to the method described in |Thirumalai| 
& Heyl (20101. Once a particular value for 7 is chosen, the 



bulk radial gas velocity at the photosphere is found using 
the relation. 



Uo = UA 



Ta 
To 



1/(7-1) 



(5) 



Thus, with these parameters now identified, a solution of 
Eq. ([l| is obtained that is continuous through the critical 
points with a tolerance for testing continuity of about 10~^, 
about an order of magnitude greater than the integration 
step size. The polytropic exponent is varied until a suitable 
set of values for the parameters {rs,Us,rf ,uf ,uo} , are ob- 
tained so that the solution would be first continuous through 
all the critical points and second, result in a temperature of 
about 700 K in the vicinity of 307?o- 



4 RESULTS AND DISCUSSION 

Observations of the circumstellar atmosphere of Betelgeuse 
indicate that dust primarily exists in the form of a shell 
at a distance of about ~ 30f?o (e.g. Bloemhof et al.||l984 |. 
However, recent observations as discussed above, seem to in- 
dicate that there may be small amounts of alumina (AI2O3) 
present rather close to the photosphere (see Perrin et al. 



2007 Verhoelst et al. 20061, within 1.57?o. However, it is 



thought that this alumina may be either transient, perhaps 
even destroyed at around 1.5J?o or even further out in the 
chromosphere (see Verhoelst et al. 20061, or is present in 



such small quantities as to be unable to support a dust- 
driven wind at these distances. There is a third possibility, 
however remote, that the small amount of alumina present 
does result in a mild dust-driven wind, but the alumina is 
transparent until it accumulates silicates on its surface (see 
Onaka et al.|1989 1 , which occurs in the dust shell at around 
307?o- 

Regardless, observations do not indicate a significant 
presence of dust in the region 1.5-Ro ~ r ;£ 30-Ro- In this pa- 
per, we address each of these possibilities, within the frame- 
work of the hybrid-dust-driven wind model. Physically, there 
are four distinct scenarios that emerge and these are listed in 
Table 2. Each of these scenarios are explored in the following 
discussion. 



Table 2. Different dust formation scenarios in Betelgeuse 

Scenario la Silicate dust forms at r = 30-Ro in reasonable 
quantities (Fd = 0.5, (5) - 1/2000) 

Scenario lb Silicate dust forms at r = 30-Ro in large 
quantities (Fd = 5, (5) - 1/200) 

Scenario 2a Alumina dust forms at r < 1.5-Ro in small 
quantities (F^ = 0.05, (5) ~ 1/20000) and 
provides nucleation sites for silicate dust 
around r = 30-Ro 

Scenario 2b Alumina dust forms at r < 1.5-Ro in large 
quantities (Fa = 0.5, (S) ~ 1/2000) and 
provides nucleation sites for silicate dust 
around r = 30-Ro 

Scenario 3a Alumina dust forms at r < 1.5-Ro in small 
quantities (F^ = 0.05, (5} ~ 1/20000) and 
is subsequently destroyed at r 7.53-Ro due 
to a chromospheric component in the atmosphere 
and silicate dust later forms at r = 30-Ro 

Scenario 3b Alumina dust forms at r < 1.5-Ro iu large 
quantities (F<j = 0.5, (<5) ~ 1/2000) and 
is subsequently destroyed at r 7.53-Ro due 
to a chromospheric component in the atmosphere 
and silicate dust later forms at r = 30-Ro 

Scenario 4a Alumina dust forms at r < 1.5-Ro in small 
quantities (F^ = 0.05, (5) ~ 1/20000) and 
is subsequently destroyed at r Ri 1.5-Ro and 
silicate dust forms at r = 30-Ro 

Scenario 4b Alumina dust forms at r < 1.5-Ro iu large 
quantities (F^ = 0.5, ((5) ~ 1/2000) and 
is subsequently destroyed at r 1.5-Ro and 
silicate dust forms at r = 30-Ro 



4.1 Scenario 1: Silicate dust forms at 30-Ro 

This is the simplest and perhaps the most likely scenario 
to fit the observations of Betelgeuse's atmosphere. Herein, 
there is no dust formation at close distances and silicate dust 
condenses at a distance of about 30-Ro where the tempera- 
ture drops to about 700 K. Stellar material is transported 
from the photosphere to this distance of 30-Ro by means of 
a Weber-Davis magneto-rotational wind. This is shown in 
Figure [1] with two cases corresponding to silicate dust for- 
mation in reasonable and large amounts at a distance of 
about 30-Ro. The gas velocities of the critical hybrid- wind 
solutions of Eq. ([T]) are shown in Figure [l] as the red and 
green long-dashed lines. These solutions start at the base 
of the wind subsonic and accelerate first through the sonic 
point at about 5.27-Ro; after this the wind is supersonic. It 
can be readily seen that the Mach numbers are small for 
the critical solution close to the photosphere (r ;S 5-Ro). Be- 
yond the sonic point, the wind mildly accelerates through 
the radial and fast Alfven points, that are nearly coincident 
upon one another, and emerges super-Alfvenic at large dis- 
tances (r > 25-Ro). A little further out, dust condenses out 
from the gas at a radial distance of rd — 30-Ro, shown by 
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dust, Scenario la 
dust, Scenario lb 
gas, Scenario la 
gas, Scenario lb 
Temperature, Scenario la 



40 
r/Ro 



Figure 1. A hybrid wind solutions are shown for Scenario 1 with 
parameters ua = 0.15vesc,o, ta = 25-Roi for different and 
remaining parameters as given in Table 1. The red and green 
solid lines trace the dust velocity profiles for Scenarios la and lb 
respectively, corresponding gas velocity profiles are shown by the 
red and green long-dashed lines. The decreasing blue solid line 
traces the temperature and should be interpreted using the right 
hand y-axis. 



the vertical dashed line marked "r^" . At this distance from 
the surface of Betelgeuse, the temperature has dropped be- 
low the dust condensation temperature for silicates of about 
~ 700 K; the temperature is given by the blue solid line and 
should be interpreted using the right-hand axis. The temper- 
ature profile shown in Figure [l] corresponds to Scenario la, 
wherein silicate dust forms in reasonable amounts at around 
30i?o- Thereafter the wind is a coupled MHD-dust-driven 
wind and the gas outflow rapidly approaches the terminal 
velocity. Thus, the solution is purely WD from the photot- 
sphere (r = Ro) to the dust condensation radius (r = r^) 
and thereafter it is hybrid. This is the primary salient fea- 
ture of the hybrid wind model developed for Betelgeuse. 

Also shown in Figure [T] are the corresponding dust ve- 
locities for these two models in this scenario, shown by the 
solid lines that lie above the long-dashed ones. It can be seen 
that the dust grains are moving radially faster than the gas 
and dragging the gas along with them. The solid red line rep- 
resents the dust velocity profile for a model with parameters 
Td = 0.5 and {5} = 1/2000; Scenario la. In order to investi- 
gate the effect of changing the average dust-to-gas ratio and 
therefore the dust parameter Td, we kept all other parame- 
ters of the model fixed, in particular, the radiation pressure 
mean efficiency and the stellar luminosity, were kept con- 
stant. Now, if the average dust-to-gas ratio is increased by 
an order of magnitude then accordingly, the dust param- 
eter Td must also correspondingly increase by an order of 
magnitude. Thus, for the second model's results shown in 



Figure [T] we took F^ = 5 and (S) = 1/200; Scenario lb. 
For this latter model, the dust velocity profile is shown with 
the green solid line that lies below the red solid line. This 
range of the dust-to-gas ratio of 1/2000 S5 {5} ^ 1/200, rep- 
resents a reasonable bound for the amount of dust in the 



atmosphere of Betelgeuse (e.g. Harper et al. 20011. It was 
seen in our earlier work (c.f. Figure 6 of Thirumalai & Heyl 
20101 that changing the dust parameter shifted the loca- 



tion of the critical points. In general, increasing the value of 
^ the dust parameter F^ results in moving the location of the 
h sonic point and fast point towards the surface of the star; 

this is the case should dust formation occur inside the sonic 

:ooo 

point in the hybrid wind model. However, it was shown ear- 
lier (see Figure 9 of Thirumalai Sz Heyl|2010 and discussion 
thereof) that formation of dust beyond the fast point does 
not influence the location of the critical points. Then, Eq. |l]) 
can simply be integrated with the presence of the Heaviside 
""o function from r — rd to r = oo for a given value of F^j. 
Thus in this case, the wind has already successfully passed 
through the critical points and emerged super- Alfvenic prior 
to dust condensation. 

Beyond about r ~ = 25-Ro, the acceleration of the 
gas in the wind due to dust drag in the second model (Fd — 
5, Scenario lb) starts to decline more steeply than in the case 
of a hybrid model with a smaller value of Td ~ 0.5 (Scenario 
la). Thus, the gas in the wind in the flrst model (Fd = 
0.5, Scenario la) at this distance, is still getting accelerated, 
therefore it's terminal velocity is slightly larger and the red 
long-dashed line lies above the green long-dashed line. Thus, 
when Fd is smaller, acceleration due to radiation pressure 
continues to have an effect, out to larger distances from the 
star. 

The effect on the dust grains is a little counter-intuitive 
and can be understood by examining Eq. Q. We can re- 
write Eq. Q by replacing the dust grain number density 
with an expression employing the dust-to-gas ratio as. 



v{r) = u{r) + 



+ 4 



6)r2 J 



\ 



1/2 



V 



(6) 



Upon examining the second term under the square-root; the 
radiation pressure term, we can see that the smaller the 
value of the average dust-to-gas ratio, the larger this term 
will be and therefore the larger the value of the dust grain 
velocity, v{r). Thus, when the radiation pressure mean ef- 
ficiency and the stellar luminosity are kept constant, then 
naturally, the dust grain velocity is larger for smaller dust- 
to-gas ratios; this is the effect seen in Figure [l] For the 
calculations carried out above, we assumed that the dust 
grains were spherical and were assumed to be on average 
(see Perrin et al. 12007 1 about 0.005/im in size with a density 
of about 4 g/cm'^. 

With regard to scattering of radiation by the dust 
grains, in the current study this was assumed to be ab- 
sent, thus precluding the complications that arise upon in- 
cluding this effect. Briefly, the inclusion of isotropic scatter- 
ing would have the effect, for the simple theory described 
here, of altering the radiation pressure mean efficiency as, 
Qrp I— >■ + , where and represent the efficien- 
cies of absorption and isotropic scattering, respectively. In 
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a more rigorous model, these could be calculated for a par- 
ticular type of dust grain and used in the equations, thus 
incorporating scattering of photons by dust grains. 

Moreover, in the framework of the current theory, the 
dust grains do not possess azimuthal velocity with respect 
to the gas. That being said, Poynting-Robertson drag due to 
scattering of radiation by dust grains would inevitably decel- 
erate the grains in the azimuthal direction, thereby altering 
the momentum equations further. In reality however, it is 
to be acknowledged that scattering is probably anisotropic 
since the dust grains may well align themselves along field 
lines. Such a detailed analysis involving the complex phe- 
nomena touched upon above, while being extremely perti- 
nent and closer to a realistic picture, was considered to be 
outside the scope of the current study, where the aim is to 
portray a simple picture. 

It is implicitly assumed in our model that dust conden- 
sation occurs abruptly at a distance where the gas temper- 
ature {Tgas) falls below the dust condensation temperature 
Tdust- In addition, we assume that that the conditions are 
conducive for grain growth. In the current, rather idealised 
treatment of the dust, we are not concerned with the radia- 
tive properties of the dust. We are therefore content with the 
assumption that the temperature of the outflow is governed 
by the gas. This was considered to be reasonable considering 
that the dust-to-gas ratio employed was small (5) = 1/2000. 
It has also been shown that the temperature profile in the 
circumstellar environment of Betelgeuse is rather compli- 
cated (e.g. Rodgers & Glassgold 19911, with effects such 



as dust-grain drag contributing to heating. Such an anal- 
ysis as well inclusion of rigorous dust radiative properties 
was considered to be outside the scope of the current paper. 
We simply assume for determining the dust-grain tempera- 
ture, that the dust grains must be in radiative equilibrium 
with the stellar luminosity field and that it is optically thin. 
Therefore, we can assume a relationship for the temperature 
profile for the dust, outside the dust condensation radius as 
(e.g. Lamers fc Cassinelli||1999 Onaka at al.|1989" l. 



Tdir) = CTo 



,V r > 307?o, 



(7) 



where, the coefficient C, depends upon the radiation pres- 
sure mean efficiency. In the current study we are treating 
the dust in a rather simplistic manner, therefore instead of 
calculating the radiation pressure mean efficiency we can as- 
sume that the coefficient C lies in a range between zero and 
unity, i.e., < C ^ 1. This then allows us to calculate a 
range of values for the dust temperature at the condensa- 
tion radius, Ttj(30-Ro). The mean of these calculated values, 
for the given range in C and the adopted values of Ro and 
To given in Table 1, is found to be {Td}\soRo ~ ^^^^ 
rather rudimentary estimate is in reasonably good agree- 



ment with the value obtained by Harper et al. (20011, of 



360 K at around 33-Ro- The dust in the current model may 
therefore be able to mimic, in a qualitative sense, the radio 



and IR behaviour of the dust model of Harper et al. (2001 1 



Given the idealistic nature of the current hybrid-MHD-dust- 
driven wind theory, the formulation of a detailed model, such 



as that described by Harper et al. (2001 1, was considered to 



be outside the scope of the current investigation. We shall 
end with a cautionary note, that a detailed analysis would 
nevertheless be required, before conclusions can be drawn 



about the dust related physics and we acknowledge this as 
a limitation of the current theory. 

In summary, the result in Figure [T] demonstrates that it 
is possible to obtain a coupling between magneto-centrifugal 
effects and the usual dust-driving mechanism. There are 
however a few pertinent observations with respect to the 
hybrid wind model in specific, that warrant mentioning. 
First, it can be clearly seen, that magneto-centrifugal ef- 
fects, with a small magnet ic field (^ 1 G ) and slow rotation, 
SI 1.2 X 10"*rad/s, (see Uitenbroek et al.^1998 ) can quite 



clearly lift material from the photosphere of Betelgeuse up 
into the circumstellar atmosphere. In the latter region, dust 
condenses in the gas to result in a hybrid wind at large dis- 
tances. This reveals that magneto-centrifugal driving can be 
an additional mechanism for solving the mystery of how to 
lift stellar material out into the circumstellar envelope in 
not only Betelgeuse, but potentially in all cool, evolved su- 
pergiant stars. Second, the temperature profile obtained in 
the solution, is reasonably consistent with observations (see 
|Plez fc Lanibert||2002{ |Lim et al.||l998| [Harper et aL||2009l ) . 
For example, the profile indicates that around 2J?o, the gas 
temperature drops to about 2840 K, lying well within the 
measured range (see Lim et al.|1998 1. 

In addition, the temperature range in the region Rq ^ 
r 5C 4.1i?o can be inferred from Figure [l]to be 1200 < T s; 
3650 K and it can be seen that lower limit for the observed 
[Fe II] emission of 2110 K (see |Harper et al. 12009] ) lies in this 
region. It can also be seen that the gas velocities are small 
in this region, consistent with the lack of observations for 
Doppler blue-shifted wind signatures. 

Third, the resulting terminal velocity is about 14 — 15 
km/s; this is consistent with present estimates (e.g. Falceta- 



Gongalves & Jatenco-Pereira 2002[ [Jura 1984) and in rea- 
sonably good agreement with the adopted value of about 
~ 10 km/s for modelling by Harper et al.| ( [2C)09| |. It is to be 
remembered that measured or inferred values for the outflow 
velocity would indicate a somewhat averaged value for gas 
velocity at a given distance from the star; however, within 
the framework of the current model, outflow velocity is cal- 
culated in the equatorial plane alone. Simulations of large- 
scale MHD convection in Betelgeuse have suggested that it 
is possible to have variation in the radial velocities on the 
order of 1 — 10 km/s in both up- and down-flowing regions 
(e.g. Dorch 2004 1 . In addition, it is also to be mentioned that 



the stellar efflux is thought to be variable in Betelgeuse, ob- 
servations seem to suggest variability over a 40 year period 
or so (see Danchi ct al. 1994), however the current frame- 
work assumes a steady-state case, thus the interpretation of 
the calculated values of outfiow obtained in this study, need 
to be tempered by these observations. In light of all this, we 
considered the agreement of velocities to be good. 

In addition, once the radial velocity profile has been 
determined, it is possible to determine the azimuthal com- 
ponent of the magnetic field. With this it becomes possible 
to estimate, a posteriori, the importance of the Lorentz force 
and any ion- neutral drag. The rationale being to check the 
validity of the ideal-MHD assumption. In the current frame- 
work, it is assumed implicitly that the ions and neutrals are 
well coupled in the equatorial plane and non-ideal MHD 
effects such as ambipolar diffusion are negligible and thus 
ignored. A way to assess this is to calculate the ratio of the 
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gyrofrequency and the momentum exchange rate (e.g. Bai 
20TTI) as, 



rrijC p(r){av)j 



(8) 



Where f3j is the ratio for the j ionic species, Zj is the 
atomic number, rrij is particle mass of the given ionic species, 
B is the magnitude of the magnetic field, e is the electronic 
charge, c is the speed of light and m is the averaged particle 
mass of the neutrals. Here p is the gas density, determined 
once the solution to Eq. ([T]) is calculated. Finally, {a-v)j is 
the ion-neutral collision momentum transfer rate coefficient 



(see Table 2.1 and Equation 2.34 of Draine 20111 and is 
given bya 

/ \ 1/2 

(av), = 2.0 X 10"' — cm' s-\ (9) 



where, /i — mjfin/{mj+pn) is the reduced mass in a typical 
ion-neutral collision, with fi„ the mean molecular mass of 
neutrals. Similarly, the electron-neutral collision momentum 



transfer rate can be determined as (e.g. Draine et al.||1983 
Bai|2011| ), 



8.3x10' 



T{r) 
100 K 



1/2 



3 - 

cm s 



'(10) 



Assuming that the largest single ionic species is Si, we can 
then calculate the ratio of the gyrofrequency to the momen- 
tum exchange rate, for both the ions (/3sj) and the electrons 

(/5e). 

We find that /Sg^" ~ 7.9 x 10"^ while ^J'^^" = 2.7 x 



lo- 



in general, we find over the entire domain for r in the 



equatorial plane, that /3si(r) ^ Pe{r) ^ 1. As mentioned 
earlier, one the central assumptions in the model is that the 
Lorentz force vanishes in the fluid, i.e., we have force- free 
MHD, therefore it is to be expected that the ratios /3 of the 
gyrofrequencies to the momentum transfer rate, would also 
accordingly be small. This then implies, that for the radial 
range considered in this study, i.e. 5^ r ^ 80iio, the Hall 
and ambipolar diffusion terms arising from electron-ion drift 
and ion-neutral drift respectively, are negligible. Addition- 
ally, the conductivity of the plasma can be estimated using 



Vish- 



the usual relati on, a = lO'^Te /^K'^/^f^-'cm"' (e 

niac et al.|2003 1, where Te is the electron temperature, such 
that Te{r) > Tgas(T'). Figure [l] shows that the minimum 
temperature for the gas far away from the photosphere is 
about Tgas" ~ 640 K. We can therefore find an estimate 
for the lower limit of the conductivity using this tempera- 
ture and estimate a value of (Jmin ~ 1.6 x 10^^ ft~^ cm"'^, 
yielding a magnetic Reynolds number (e.g. Davidsonl|2001 1 
of Rem = ULpoCmin ~ 10^', SO Ohmic diffusion is unim- 
portant {U ~ 10 km/s and L « 7?o are typical velocity and 
length scales for Betelgeuse and po is the permeability of the 
vacuum, Atv x 10~^f2s/m). The assumption of effectively in- 
finite plasma conductance, is then a reasonable assumption. 
Thus, we see that the central assumption of ideal-MHD in 
the equatorial plane of Betelgeuse is a reasonable first ap- 
proximation to make, for the purpose of conveying a simple 
picture. 

Finally, the Mach numbers in the inner region of the 
circumstellar envelope are small and actually, do not ex- 
ceed unity before the sonic point, located at around 5.27i?o. 
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Figure 2. The azimuthal velocity of the gas is shown as a function 
of the radial distance (red solid line) for Scenario la. The inset 
shows a magnified region around the peak of the profile around 
the sonic point. 



Therefore our model does not suffer from the same draw- 
backs as many Alfven wave models that have the artifact of 
having large Mach numbers close to the surface of the star. 
Indeed, in our model the wind only becomes super-Alfvenic 
beyond the Alfven point at around 25-Ro and even then, only 
mildly so. 

In Figure[2]we have shown the azimuthal velocity, u^{r) 
of the gas (red solid line). The dust is assumed to co-rotate 
with the gas, thus the azimuthal velocity profile shown in 
the figure is also true for the dust (for r ^ r^j), since there 
is no drag between the dust and the gas in the azimuthal 
direction. The azimuthal velocity profile is obtained once 
the radial velocity profile is known (see [Thirumalai fc Heyl| 



2010 1. The velocity profile obtained is typical for a magneto- 
rotational wind. The dotted lines in the figure and in its in- 
set represent the rotation of the star (see [Uitenbroek et ah] 
1998). It is readily seen that the azimuthal wind velocity 
closely traces the observed rotation of the star and only be- 
gins to depart markedly after r ~ 4.5i?o. This decoupling 
from stellar rotation occurs in the inner wind region of the 
magneto-rotational wind, at a physical distance correspond- 
ing to « 2.03 X lO'^^cm, for Betelgeuse. It is also interest- 
ing to note that at large distances the azimuthal velocity of 
the gas is strongly de-coupled from the stellar rotation rate. 
Therefore, the chromospheric component in the atmosphere 
at a few stellar radii may also appear to be de-coupled from 
stellar rotation as was observed by ( Harper fc Brown|2006 |. 
However, it is to be remembered that the hybrid-MHD-dust- 
driven theory presented in this study, concerns itself with the 
equatorial plane alone, while the observations of Harper fc| 



Brown ( 2006 1 are spatially resolved in the cross-dispersion 



direction; this may have an averaging effect along the dis- 
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To = 3650 K 
7], = 3000 K 
T„ = 2600 K 



MARCS model 
Vcrhoclst et al (2006) 



MOLspherc model, Perriu et al (2007) 



2.0 2.5 

r/Ro 



Figure 3. Temperature profile in tlie close circumstellar environ- 
ment with two photosplieric spots of temperatures 2600 K and 
3000 K. Tile two dasiied liorizontal lines sliow calculated alumina 
condenstation temperatures from two different models. The ver- 
tical green shaded region shows the MOLsphere region predicted 
by Perrin et al (2007). Only the model with photospheric spot 
temperature of Tq = 2600 K has an overlap with the MOLsphere 
region for the alumina condensation temperature range. 



persion direction. Thus, the theoretically obtained rotational 
decoupling that occurs in the equatorial plane, in the current 
study, may well be quite a bit different from the observations 



of Harper & Brown ( 2006 1 . Therefore any inferences drawn 



in this regard must be tempered by the observation that the 
current model is a rather simple picture, with concomitant 
limitations. 

We now turn our attention to the question of spots on 
the surface of Betelgeuse and the related question of tem- 
perature inhomogeneities. Recent observations (see |Wilson| 



et al.||1992[ |1997| [Buscher et al.|p90l [Haubois et ai.||2009 1 

suggest that there may be a few (on the order of 2 or 3) large 
spots on the surface of Betelgeuse, indicating that there are 
temperature inhomogeneities on the surface of Betelgeuse. 
We therefore chose to model Betelguese with two spots on 
its equator (see |Soker k, Clayton|1999l ). 

In the current study however, these spots were taken to 
be colder than the effective temperature so that the tem- 
perature profile close to the photosphere 1.0 55 r ^ 3.5-Ro 
may be investigated so as place constraints on the alumina 
formation region, if at all, close to the star. The spots were 
considered to have temperatures of 2600 K and 3000 K re- 
spectively, well within fluctuations of about 1000 K or so, 



about the effective temperature (see Schwarzschild 1975 1 



The magnetic field was left unchanged so that the effect of 
changing a single parameter could be patently established. 
The solitary requisite constraint that was placed however, 
was that regardless of the photospheric spot temperature, 



the temperature at a distance of about 30-Ro should be in the 
vicinity of 700 K, for silicate dust formation in the circum- 
stellar shell at this distance. Thus, the radial Alfven point 
was left unchanged. Therefore the other free parameter that 
was varied in order to achieve an efflux was the bulk gas 
radial velocity at the surface, uq. The gas velocity profiles 
were then solved for by integrating Eq. ([T|, above the spots 
located at the photosphere. Figure |3] shows the temperature 
profiles obtained in the wind ahead of the spots in the close 
circumstellar environment of Betelgeuse. The temperature 
profile obtained using the effective temperature of 3650 K is 
also shown for comparison as the top most blue solid line. 
It can be seen that the lower the photospheric temperature, 
the flatter the temperature profile is, in the inner wind re- 
gion. Thus the red solid line (Tq = 3000 K) lies below the 
blue solid line (To = 3650 K) and above the green solid line 
(To = 2600 K). 

Also sho wn in Figure [3| are two temperature bounds 
calculated by Verhoelst et al. ( 2006 1 for alumina dust con- 
densation. The upper bound of 2400 K is the value that 
they arrived at by demanding radiative equilibrium of dust 
grains using the radiative transfer code MODUST, for alu- 
mina condensation. The lower bound of 1900 K on the 
other hand, is the temperature that their MARCS model 
predicts in the region about 0.5i?o above the photosphere. 
It is to be kept in mind however, that at such short dis- 
tances, the MARCS model may have limitations due to 
acoustic wave heating and the breakdown of the assump- 
tion of local thermal equilibrium. With this caveat in place, 
it can be seen that this range of temperature only occurs 
in the region 2.55i?o ^ ^ 3.16i?o for the model with 
To = 3650 K, indicating the possible alumina dust con- 
densation region. Similarly, the corresponding range for the 
model with the spot of temperature To = 3000-/? places this 
region to be 1.99i?o ^ ^ 2.83-Ro. Both these regions lie 
well outside the observed thin shell for the MOLsphere of 
1.33i?o ^ r- < 1.5i?o (e.g. [Perrin et al.||2007| [Verhoelst et al. 
|2006| ; shown as the green band in Figure |3j 

It is however interesting to note, that the model for the 
photospheric spot temperature of 2600 K, predicts a possible 
alumina condensation region of 1.37iio ^ 5; 2.37-Ro- This 
region quite interestingly, overlaps slightly with the MOL- 
sphere model of Per rin et al.^ ( ,2007_) . However, it is to be re- 
membered that current observations do not reveal any dust 
in the region 1.57?o ^ r ^ 20 — 307?o (e.g. Verhoelst et al. 
[2006] ). Thus, we see that within the framework of the hybrid- 
MHD-dust-driven model presented in this paper for Betel- 
geuse, it is possible to form alumina dust very close to the 
photosphere in the region of interest; 1.33i?o ^ r 1.5i?o, 
by lowering the photospheric temperature to about 2600 K. 
For temperatures above this, our model predicts the possi- 
ble alumina condensation region to lie further out from the 
photosphere and there is no overlap between our model and 
the MOLsphere model of [Perrin et al.| ( |2007[ ). 

Additionally, it can be seen in Figure [3j that tempera- 
tures predicted by the models with photospheric spots are 
different at the same height above the photosphere, in the 
close circumstellar environment of Betelgeuse. Thus, the ob- 
served variability in temperature (see Lim et al.|1998 1 at the 
same radial distance from the photosphere can, at least in 
part, be attributed to the presence of spots. The temper- 
ature profile in the wind above the spots is different from 
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that above normal regions. Thus, our model with equato- 
rial spots, can reproduce these differences in temperature, 
at least qualitatively, indicating an asymmetric temperature 
distribution. 

Presently, we discuss the other possible scenarios for 
dust formation listed in Table 2. 

The motivation for these scenarios is the observation 
that alumina dust may be present in small quantities close 
to the photosphere of Betelgeuse. Therefore, one can nat- 
urally ask the question whether alumina first forms in suf- 
ficient amounts to facilitate a slow dust-driven wind and 
second, what would happen if alumina dust is subsequently 
destroyed at some distance from the star? In this regard 
there are a few possibilities. It is to be kept in mind that 
there are no observations that support the presence of dust 
between ~ 1.5-Ro and about 20 — 30-Ro (e.g. Verhoelst et al 



20061. Therefore, the first possibility is we can assume that 



alumina forms in both small and sufficient quantities close 
to the star {Scenanos 2a and 2h respectively), to facilitate 
a slow wind, but the alumina dust is transparent until it 
accumulates silicates on its surface beyond 30i?o. In both 
these models, the assumption is that alumina dust provides 
nucleation sites for silicates to form upon. 

The second possibility is that alumina dust forms close 
to the star (in small and large quantities; explored in Sce- 
narios 3a and 3b respectively) . This results in a slow MHD- 
dust-driven hybrid wind. The dust is then destroyed at 
around « 7.53-Ro due to a chromospheric component in the 



atmosphere where the temperatures are high (e.g. Lobel & 
|Dupree|[200T| . Note that in ( |Lobel fc Dupree|[2001| ) they 
adopted a value for the radius of the photosphere as 7OO-R0 , 
whereas we have adopted 650/?© , hence we obtain the outer 
limit for the chromosphere as ~ 7.53-Ro rather than 7Ro- 
Silicate dust later condenses at large distances ~ 30-Ro- It is 
implicitly assumed in this model that alumina dust is trans- 
parent so as not to reveal any dust signature between 1.5-Ro 
and 30-Ro- 

The third possibility is we can assume formation of alu- 
mina in the region r ^ 1.5-Ro (in small and large quantities; 
explored in Scenarios 4<i and 4b respectively). This dust 
then gets destroyed at a distance of 1.5-Ro due to perhaps 
convective turbulence in regions closer to the photosphere 
and changes in pressure or perhaps due to temperature vari- 
ability in the chromosphere. Whatever the reason may be, 
dust is not seen between 1.5-Ro and about 20 — 307?o. Silicate 
dust later condenses at large distances. 

The purpose of these scenarios is not to elaborate on the 
details of dust spallation, but rather to ask the pertinent 
question that, within the framework of the hybrid-MHD- 
dust-driven wind theory, is it possible to achieve an efflux, 
should alumina first form and perhaps even be destroyed 
in the wind, at some distance? In addition, it is to be kept 
in mind that for the models presented in Scenarios 2 — 4 
the photospheric temperature was considered to be To — 
2600 K. 



4.2 Scenarios 2 and 3: Alumina forms and has an 
influence on the wind 

Figure |4] shows both Scenarios 2 and 3 each containing two 
hybrid wind models. Scenario 2 is shown in the upper panel 
of Figure |4] while Scenario 3 is shown in the lower panel. The 




Figure 4. Dust and gas velocity profiles for Scenarios 2 (upper 
panel) and 3 (lower panel). The dust velocities are shown using 
solid lines and the gas velocities with long-dashes lines. The dust 
and gas velocity profiles for the models with greater dust forma- 
tion, i.e. larger F^, lie above those with lesser amount of alumina 
dust formation in the inner wind region. Scenario 2 shows infiu- 
ence of the presence of both alumina and silicate dust species in 
the wind, whereas, Scenario 3 shows the effect of alumina spalla- 
tion at the edge of the chromosphere at 7.53-Ro- 



red and green solid lines represent the dust velocity profiles 
in these scenarios for small and large amounts of alumina 
condensation, respectively. Similarly, the long-dashed green 
and red lines represent the corresponding gas velocity pro- 
files. Quite importantly, in Scenarios 2a and 2b, the numer- 
ator in Eq. ([2| was modified to accommodate for the two 
different dust species; alumina and silicates as, 



2\l-7 



Sg 

X 



N{w,x) = (2'yST{wx'^f 

X (1 - r:^^™'"" ■ e{x - xll'') - Ff 

x{wx'^ - 1)^ + Snx'^iw - 1) (1 - Swx"^ + {wx^ + l)w) , (11) 



where, x^J'^ — 1.38 for alumina dust condensation and a;^^' — 
30 for silicate dust condensation. Thus, Eq. ([TJ is integrated 
as a pure WD wind in the region Rq ^ r ^ 1.38-Ro, then with 
alumina dust in the region 1.38-Ro ^ r 5C 30-Ro and finally 
with both alumina and silicate dust in the region 30-Ro ^ 
r 5C 40-Ro. In Scenario 2 (upper panel of Fig. [4| the portion 
of the hybrid solution in the region 1.38-Ro ^ ^ 30-Ro is a 
critical solution as it navigates through the critical points. 
The portions of the solutions outside this interval are not 
part of the critical solution, but rather lie elsewhere in the 
u — r phase space of each of the models, respectively. 

The alumina condensation radius of 1.38-Ro lies in the 



range suggested by Perrin et al. (20071, for the thin molec- 



ular shell that contains alumina in the close circumstellar 
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environment of Betelgeuse. This value also lies inside the 
green shaded region of Figure |3] The lower red solid lines in 
Figure |4] in both the upper and lower panels, represent the 
dust velocity profiles for models with formation of alumina 
dust in small quantities {{S) = 1/20000 ). This is an order 
of magnitude less than what is expected for silicate dust in 
the dust shell at around 30iio- Accordingly, we assumed an 
order of magnitude smaller value for the dust parameter, 
Fd = 0.05; it is to be mentioned that even smaller values 
may be adopted. 

The upper green solid lines in both panels, prior to the 
silicate condensation region, represent the dust velocity pro- 
files for models with formation of alumina dust in large quan- 
tities ( (S) = 1/2000). The dust parameter was then given by 
Td = 0.5. The purpose being to explore a somewhat extreme 
case in which alumina forms in large enough amounts to 
sustain a mild MHD-dust-driven hybrid wind, but is trans- 
parent to observations. While such a scenario may well be 
remote, we explore it for the sake of completeness. 

It can clearly be seen in both Scenarios 2 and 3, as ex- 
pected, that if the dust parameter is larger, then it results in 
a greater acceleration of the wind closer to the photosphere. 
Thus, the green lines lie above the red ones. For smaller val- 
ues of the dust parameter, on the other hand, the effect of 
this acceleration due to dust, continues to have an effect out 
to larger distances. As a result, the red long-dashed lines 
rise above the green long-dashed lines at some distance. Ul- 
timately the terminal velocity of the models with the smaller 
value of Fd are slightly larger due to this effect. 

In the upper panel in Figure [4] at the silicate dust con- 
densation radius of 30i?o, we notice a discontinuity, since 
the silicate dust condenses as well. The red solid line in 
this region corresponds to a model with parameters Td ~ 

palumi„a_^pBilicatc ^ 0.05-^0.5 = 0.55 with (5total) = 1/2000. 

Thus Eq. |T| is integrated outwards from r — rd with these 
parameters. Similarly the green solid line in this region in 
the upper panel, represents a model with parameters corre- 
sponding to Fd = r^i"'"i"'' + rf' = 0.5 0.6 = 1.1, again 
with (5totai) = 1/2000. Notice that we have kept the dust-to- 
gas ratio fixed in both cases to the same value. The silicate 
dust condenses from the gas abruptly at 30i?o. As a result, 
the discontinuity in the dust velocity across the silicate dust 
condensation radius is seen in both models. Ifowever, the 
jump in the velocity is much larger for the green solid line 
in comparison to the red solid line as this model has a higher 
value of Fd for the same (5totai). 

On the other hand, in the lower panel of Figure [4] show- 
ing scenarios 3a and 3b, alumina dust is assumed to form 
close to the photosphere, but is then destroyed at a distance 
of about 7.53J?o due to a chromospheric component. There 
are some rather subtle differences from Scenarios 2a and 2b 
and these are detailed below. 

In the lower panel of Figure |4] once again the lower 
solid red line shows the dust velocity profile for a model 
with parameters f^'"'"'"'' = 0.05 and ((Jaiumina) = 1/20000. 
The upper green solid line meanwhile, represents the dust 
velocity profile for a model with parameters r|^'"™'"'' = 0.5 
and {(Jaiumina) = 1/2000. As described earlier, the solutions 
start off at the base of the wind at the photosphere, subsonic 
and alumina dust formation occurs at r = 1.38-Ro. Thus 
from r = Ro to r = 1.38-Ro, Eq. ([T]| is integrated as a pure 
WD wind without the dust parameter. It is to be mentioned 



that in this region the solutions are not part of the pure 
WD critical solution. Alumina dust condensation radius is 
shown at 1.387?o for the two models in this scenario. As 
can be seen, the model with the larger dust parameter has 
a greater velocity, evidenced by the green lines lying above 
the red ones. 

In the region l.SSiio ^ r ^ 7.53-Ro the winds are hy- 
brid MHD-dust-driven winds and Eq. ([T]) is integrated with 
the presence of dust parameters for each of the two models. 
In this region the solutions for the gas velocities are parts 
of the critical solutions for each of the two hybrid-MHD- 
dust-driven winds. Alumina dust spallation is assumed to 
occur at a distance of about ~ 7.53-Ro and this produces 
the sharp discontinuity in both the dust and gas velocity 
profiles seen in the lower panel of Figure |4] In the region 
7.53-Ro S r ^ 30i?o, Eq. ([l| is therefore integrated with- 
out the dust parameter, i.e. as a pure WD wind, for the 
two models. However, the gas velocities in this latter re- 
gion, are no longer part of the respective critical solutions. 
Hence, these solutions pass through only the radial Alfven 
point, but not the fast point. Thus at large distances, these 
solutions are sub-Alfvenic. Therefore they constitute the so- 
called failed wind solutions as r — >■ oo. As a result these 
models become theoretically non-viable as one of the pri- 
mary requirements is that the wind should be super- Alfvenic 
at large distances. Thus with alumina dust spallation in be- 
tween the sonic point and the radial Alfven point, at around 
7.53-Ro, we see that it is not possible to sustain a physically 
viable efflux. 

For the sake of completeness, we have included silicate 
dust in Scenario 3. At a distance of 30iio in the wind, silicate 
dust condensation occurs and beyond 30iio, the wind is once 
more a hybrid MHD-dust-driven wind. The dust velocities 
are once again calculated after the gas velocity profiles are 
determined, according to Eq. Q. It can be seen that the 
gas outflow velocity profiles for the two models cross one 
another at around the radial Alfven point. As a result the 
dust and gas velocity profiles in the region r ^ 30-Ro, for 
the model with greater amount of alumina (green dashed 
and solid lines. Scenario 3b) lie slightly below the one with 
the lesser amount of alumina (red dashed and solid lines. 
Scenario 3a). In both these models it is assumed that silicate 
dust forms with a dust-to-gas ratio of (^silicate) = 1/2000, 
with Fd = 0.5. 

In summary we see in the upper panel of Figure |4j that 
it is possible to obtain a hybrid-MHD-dust-driven wind if 
alumina forms close to the photosphere and results in a mild 
wind transporting stellar material to large distances. Silicate 
dust condenses at a distance of 30-Ro and adds to the stellar 
efflux. Inspecting the upper panel of Figure |4] reveals that 
the model with a larger amount of alumina formation close 
to the star (Scenario 2b), is not viable, since observations do 
not indicate such a large amount of alumina close to the pho- 
tosphere. In the second model (Scenario 2a) with a smaller 
value of Fd, while the amount of dust may well be reason- 
able, current observations do not reveal any dust velocities. 
The only way we can reconcile this model with observations 
is by assuming that the alumina dust is transparent. 

In the lower panel of Figure |4] on the other hand, once 
again it is seen that it is possible to have alumina formation 
at close distance, within 1.5-Ro and subsequent spallation at 
around ~ 7.53-Ro and still achieve a solution with the forma- 
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tion of silicate dust. However, as stated earlier, the gas veloc- 
ity solutions are not part of the critical solution. As a result, 
these solutions are sub-Alfvenic as r — > oo. Therefore, they 
cannot be considered to be viable wind solutions. As a result, 
alumina spallation at a distance of 7.53i?o does not result 
in a wind, within the framework of the current steady-state 
hybrid-MHD-dust-driven theory. In addition, the spallation 
of alumina dust at a distance of around ~ 7.537?o results in 
the discontinuity seen in the lower panel of Figure [4] Such 
large discontinuities are a source of concern, particularly be- 
cause they can lead to fluid flow instabilities and shocks. 
However, this is not within the realm of investigation of the 
current theory, a simple steady-state treatment. It may be 
possible that such instabilities may precipitate in motion of 



material back towards the photosphere as is seen in ( Lobel 
fc Dupree|[200l] ). Thus we see that scenarios 3a and 3b are 
theoretically possible; however, the primary concern is that 
they predict sub-Alfvenic velocities at large distances. Thus 
the results of Scenarios 3a and 3b cannot be completely rec- 
onciled with observations. This may well point towards the 
inference that the simplest scenario discussed earlier. Sce- 
nario la, may well the be the most pertinent. 

We now turn our attention to the final scenario before 
summarising the results obtained in this study. 



4.3 Scenario 4: Alumina spallation close to the 
photosphere 

Presently, we discuss the cases where alumina dust conden- 
sation occurs in both small and large quantities close to 
the photosphere and subsequent spallation occurs at around 
1.5-Ro; Scenarios 4a and 4b respectively. The rationale being, 
that the reason why alumina dust is not seen between 1.5-Ro 
and around 30i?o is that it gets destroyed at around 1.5-Ro- 
This is an alternative to assuming that the alumina grains 
are transparent until they accumulate silicates on their sur- 
faces. The question then is whether it is possible to achieve a 
hybrid outflow with alumina dust formation and spallation 
included in the framework of the hybrid MHD-dust-driven 
picture. It is to be mentioned at the very outset that the 
formation of alumina in both small and large quantities was 
not seen to produce any appreciable influence on the stellar 
efflux in this scenario. The silicate dust and corresponding 
gas velocity profiles obtained were seen to be nearly identical 
to those shown in Figurejl] Thus, formation and subsequent 
spallation of alumina, all within r = 1.5i?o was not seen 
to be a physically relevant picture for the stellar efflux of 
Betelgeuse. Nevertheless, a brief description of the calcula- 
tion and some pertinent points are conveyed below. 

For Scenarios 4a and 4b it was assumed that alumina 
condensation occurs at r = 1.38-Ro in small and large quan- 
tities. Thus Scenario 4a has parameters r;^'"™'"'" = 0.05 
and {(5aiumina) = 1/20000, while Scenario 4b has parame- 



ters r^'™'"^ = 0.5 and {(^alumina) = 1/2000. From r = Ro to 
r = 1.38iio, Eq. |T| is integrated as a pure WD wind with- 
out the dust parameter. It is to be mentioned that these 
parts of the solutions are not part of the WD critical solu- 
tion. In the region 1.38i?o ^ r ^ 1.5-Ro, the winds are hybrid 
MHD-dust-driven winds and Eq. ([T]) is integrated with the 
presence of dust parameters for the two models. However 
alumina is dust is then abruptly destroyed at a distance of 
1.5-Ro and therefore in the region 1.5-Ro < ?" < 30-Ro the 



wind is a pure WD type of wind and Eq. |T| is integrated 
without the dust parameter. The portion of the solution in 
the region 1.5-Ro < r < 30-Ro corresponds to a pure WD 
critical solution, one that passes through all three critical 
points. The critical solution then carries the wind out to a 
distance of 30-Ro at which point, silicate condensation occurs 
and silicate dust grains form. The dust and gas velocities for 
the two models in Scenario 4a and 4b in the region r 5C 1.5-Ro 
were calculated to be miniscule; ^ 10~^° cm/s and as a re- 
sult immeasurable and not detectable in the form of wind 
signatures. Thus in this scenario, for all intents and pur- 
poses, alumina dust condensation can be considered to be 
absent in the wind as it does not produce any appreciable 
effect, in which case Scenario 4 quite simply reduces to Sce- 
nario 1; the most viable picture thus far in the discussion. 

This brings us to the end of our discussion of the results 
obtained in this study. In the following section, the findings 
are summarised. 



5 CONCLUSION 

In this paper we presented a hybrid-MHD-dust-driven wind 
model for the red supergiant Betelgeuse (a— Orionis). The 
model is a direct application of our previously derived the- 
ory that consists of incorporating a dust-driven wind with a 



Wober-Davis MHD equatorial wind (see Thirumalai & Heyl 
[2010) . 

Overall, the results shown above indicate that MHD 
effects and radiation pressure on dust grains can have a 
complementary role to play in the winds of supergiants 
such as Betelgeuse, alongside the altogether complicated and 
involved physics of thermal pulsation and convection and 
other equally complex phenomena such as MHD or acoustic 
waves. 

Within the framework of the model, we investigated 
four different scenarios for dust formation in the atmosphere 
of Betelgeuse. It was seen that the simplest hybrid wind sce- 
nario; Scenario la, was perhaps the most viable one as well. 
In this picture, a pure WD wind was assumed to begin at 
the surface of the star, one that would eventually leave the 
star as a hybrid-MHD-dust-driven wind after the formation 
of dust grains at the dust condensation radius at around 
30-Ro. This provided a mechanism for lifting stellar material 
from the photosphere of the star up into the atmosphere 
while maintaining low Mach numbers for the wind in the in- 
ner wind region. The dust formation was assumed to occur 
abruptly at the dust condensation radius where the tem- 
perature is low enough for silicate grains to condense. In 
our model, all the dust grains were assumed to be perfectly 
spherical with identical size. It was implicitly assumed that 
radiation pressure was purely in the radial direction without 
scattering. The opacity of the grains were implicitly assumed 
to be such that all of the radiation impinging on the grains 
was absorbed and imparted momentum to the grains. The 
resulting drag force was assumed to be purely radial as well. 

In Scenario la, since dust grain condensation occurs 



outside the fast point, as in our earlier work (see Thirumalai 
fc Heyl|2010 1, it was seen that adopting a range of value for 
values the dust parameter Fd, resulted in different hybrid 
winds with different terminal velocities. 

It was also seen that by adopting the values for different 
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parameters for Betelgeuse, as given in Table 1, the result- 
ing hybrid wind model is able to rationalise several of the 
observed features of this red supergiant, such as the pre- 
dicted wind velocities, observed atmospheric temperatures 
at different distances from the photosphere, observed region 
of dust formation and temperature inhomogeneities in the 
circumstellar envelope. 

We additionally investigated the possibility of having 
more than one dust species form in the wind. Scenarios 2-4 
dealt with the formation of alumina dust close to the photo- 
sphere of Betelgeuse. We also investigated the effect on the 
efflux, should alumina be destroyed at some distance. In each 
of these scenarios it was seen that it was possible to form 
a hybrid-MHD-dust-driven wind, however it was not possi- 
ble to reconcile the predictions of these models completely 
with the current observations regarding dust formation and 
lack of wind signatures interior to the dust shell of Betel- 
geuse. In addition, we found that for these scenarios to work, 
within the framework of the hybrid-MHD-dust-driven wind 
theory, it was required that the surface temperature should 
be around To = 2600 K, to overlap with recent observations 
for the close circumstellar environment of Betelgeuse (e.g. 
Perrin et al. 20071. Therefore, it is to be concluded that 



while it may well be possible for alumina to form close to 
the photosphere of Betelgeuse, it probably does not influ- 
ence the wind outflow interior to ~ 30i?o- However, it is to 
be remembered that the current model is merely a steady- 
state treatment of an infinitely more complicated picture. 
All that can be said at this stage is that the current theory 
seems to point towards the simplest case; Scenario la, as the 
most likely candidate. It is to be mentioned that the current 
model may play a complementary role alongside the alto- 
gether different and complex mechanisms for stellar outflow 
involving MHD or acoustic waves. 

Future work will extend the modelling towards a more 
realistic full two- or three-dimensional picture, concomi- 
tantly including a poloidal magnetic field and allowing the 
entire system to be described dynamically with time varying 
magnetic and velocity fields with convection. The final and 
indeed ultimate picture would of course, be to then include 
the stochastic effects of dust formation and dynamics, along- 
side calculation of the radiative properties of the dust, as is 
done in some recent work (see Woitke Sz Niccolini|2005 Lo- 
bel fc Blomme|2008 1. Realisation of this task will ultimately 
explore the true nature of the outflow from Betelgeuse and 
other red supergiants, which is now beyond doubt, under- 
stood to be more complicated than at first imagined and 
yet, presenting a hurdle surmountable in steps. Perhaps the 
very first of those steps, is the steady-state model described 
in the current study. 
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